SARS-CoV-2 mutations: the biological trackway towards viral fitness.

The outbreak of pneumonia-like respiratory disorder at China and its rapid transmission world-wide resulted in public health emergency, which brought lineage B betacoronaviridae SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) into spotlight. The fairly high mutation rate, frequent recombination and interspecies transmission in betacoronaviridae are largely responsible for their temporal changes in infectivity and virulence. Investigation of global SARS-CoV-2 genotypes revealed considerable mutations in structural, non-structural, accessory proteins as well as untranslated regions. Among the various types of mutations, single-nucleotide substitutions are the predominant ones. In addition, insertion, deletion and frame-shift mutations are also reported, albeit at a lower frequency. Among the structural proteins, spike glycoprotein and nucleocapsid phosphoprotein accumulated a larger number of mutations whereas envelope and membrane proteins are mostly conserved. Spike protein and RNA-dependent RNA polymerase variants, D614G and P323L in combination became dominant world-wide. Divergent genetic variants created serious challenge towards the development of therapeutics and vaccines. This review will consolidate mutations in different SARS-CoV-2 proteins and their implications on viral fitness.

In SARS-CoV-2 delta variants, Spike-P681R and D950N promote membrane fusion, Spike-P681R enhances spike cleavage, but neither substitution affects pathogenicity in hamsters.

BACKGROUND: The SARS-CoV-2 delta (B.1.617.2 lineage) variant was first identified at the end of 2020 and possessed two unique amino acid substitutions in its spike protein: S-P681R, at the S1/S2 cleavage site, and S-D950N, in the HR1 of the S2 subunit. However, the roles of these substitutions in virus phenotypes have not been fully characterized. METHODS: We used reverse genetics to generate Wuhan-D614G viruses with these substitutions and delta viruses lacking these substitutions and explored how these changes affected their viral characteristics in vitro and in vivo. FINDINGS: S-P681R enhanced spike cleavage and membrane fusion, whereas S-D950N slightly promoted membrane fusion. Although S-681R reduced the virus replicative ability especially in VeroE6 cells, neither substitution affected virus replication in Calu-3 cells and hamsters. The pathogenicity of all recombinant viruses tested in hamsters was slightly but not significantly affected. INTERPRETATION: Our observations suggest that the S-P681R and S-D950N substitutions alone do not increase virus pathogenicity, despite of their enhancement of spike cleavage or fusogenicity. FUNDING: A full list of funding bodies that contributed to this study can be found under Acknowledgments.

T-Toxin Virulence Genes: Unconnected Dots in a Sea of Repeats.

In 1970, the Southern Corn Leaf Blight epidemic ravaged U.S. fields to great economic loss. The outbreak was caused by never-before-seen, supervirulent, Race T of the fungus Cochliobolus heterostrophus. The functional difference between Race T and O, the previously known, far less aggressive strain, is production of T-toxin, a host-selective polyketide. Supervirulence is associated with ~1 Mb of Race T-specific DNA; only a fraction encodes T-toxin biosynthetic genes (Tox1). Tox1 is genetically and physically complex, with unlinked loci (Tox1A, Tox1B) genetically inseparable from breakpoints of a Race O reciprocal translocation that generated hybrid Race T chromosomes. Previously, we identified 10 genes for T-toxin biosynthesis. Unfortunately, high-depth, short-read sequencing placed these genes on four small, unconnected scaffolds surrounded by repeated A+T rich sequence, concealing context. To sort out Tox1 topology and pinpoint the hypothetical Race O translocation breakpoints corresponding to Race T-specific insertions, we undertook PacBio long-read sequencing which revealed Tox1 gene arrangement and the breakpoints. Six Tox1A genes are arranged as three small islands in a Race T-specific sea (~634 kb) of repeats. Four Tox1B genes are linked, on a large loop of Race T-specific DNA (~210 kb). The race O breakpoints are short sequences of race O-specific DNA; corresponding positions in race T are large insertions of race T-specific, A+T rich DNA, often with similarity to transposable (predominantly Gypsy) elements. Nearby, are 'Voyager Starship' elements and DUF proteins. These elements may have facilitated Tox1 integration into progenitor Race O and promoted large scale recombination resulting in race T. IMPORTANCE In 1970 a corn disease epidemic ravaged fields in the United States to great economic loss. The outbreak was caused by a never-before seen, supervirulent strain of the fungal pathogen Cochliobolus heterostrophus. This was a plant disease epidemic, however, the current COVID-19 pandemic of humans is a stark reminder that novel, highly virulent, pathogens evolve with devastating consequences, no matter what the host-animal, plant, or other organism. Long read DNA sequencing technology allowed in depth structural comparisons between the sole, previously known, much less aggressive, version of the pathogen and the supervirulent version and revealed, in meticulous detail, the structure of the unique virulence-causing DNA. These data are foundational for future analysis of mechanisms of DNA acquisition from a foreign source.

A computational approach to biological pathogenicity.

The current pandemic (COVID-19) has made evident the need to approach pathogenicity from a deeper and more systematic perspective that might lead to methodologies to quickly predict new strains of microbes that could be pathogenic to humans. Here we propose as a solution a general and principled definition of pathogenicity that can be practically implemented in operational ways in a framework for characterizing and assessing the (degree of) potential pathogenicity of a microbe to a given host (e.g., a human individual) just based on DNA biomarkers, and to the point of predicting its impact on a host a priori to a meaningful degree of accuracy. The definition is based on basic biochemistry, the Gibbs free Energy of duplex formation between oligonucleotides and some deep structural properties of DNA revealed by an approximation with certain properties. We propose two operational tests based on the nearest neighbor (NN) model of the Gibbs Energy and an approximating metric (the h-distance.) Quality assessments demonstrate that these tests predict pathogenicity with an accuracy of over 80%, and sensitivity and specificity over 90%. Other tests obtained by training machine learning models on deep features extracted from DNA sequences yield scores of 90% for accuracy, 100% for sensitivity and 80% for specificity. These results hint towards the possibility of an operational, objective, and general conceptual framework for prior identification of pathogens and their impact without the cost of death or sickness in a host (e.g., humans.) Consequently, a reasonable prediction of possible pathogens might pave the way to eventually transform the way we handle and prepare for future pandemic events and mitigate the adverse impact on human health, while reducing the number of clinical trials to obtain similar results.

Myosin-II proteins are involved in the growth, morphogenesis, and virulence of the human pathogenic fungus Mucor circinelloides.

Mucormycosis is an emerging lethal invasive fungal infection. The infection caused by fungi belonging to the order Mucorales has been reported recently as one of the most common fungal infections among COVID-19 patients. The lack of understanding of pathogens, particularly at the molecular level, is one of the reasons for the difficulties in the management of the infection. Myosin is a diverse superfamily of actin-based motor proteins that have various cellular roles. Four families of myosin motors have been found in filamentous fungi, including myosin I, II, V, and fungus-specific chitin synthase with myosin motor domains. Our previous study on Mucor circinelloides, a common pathogen of mucormycosis, showed that the Myo5 protein (ID 51513) belonging to the myosin type V family had a critical impact on the growth and virulence of this fungus. In this study, to investigate the roles of myosin II proteins in M. circinelloides, silencing phenotypes and null mutants corresponding to myosin II encoding genes, designated mcmyo2A (ID 149958) and mcmyo2B (ID 136314), respectively, were generated. Those mutant strains featured a significantly reduced growth rate and impaired sporulation in comparison with the wild-type strain. Notably, the disruption of mcmyo2A led to an almost complete lack of sporulation. Both mutant strains displayed abnormally short, septate, and inflated hyphae with the presence of yeast-like cells and an unusual accumulation of pigment-filled vesicles. In vivo virulence assays of myosin-II mutant strains performed in the invertebrate model Galleria mellonella indicated that the mcmyo2A-knockout strain was avirulent, while the pathogenesis of the mcmyo2B null mutant was unaltered despite the low growth rate and impaired sporulation. The findings provide suggestions for critical contributions of the myosin II proteins to the polarity growth, septation, morphology, pigment transportation, and pathogenesis of M. circinelloides. The findings also implicate the myosin family as a potential target for future therapy to treat mucormycosis.

SARS-CoV-2 Delta Variant: Interplay between Individual Mutations and Their Allosteric Synergy.

Since its first appearance in April 2021, B.1.617.2, also termed variant Delta, catalyzed one major worldwide wave dominating the second year of coronavirus disease 2019 (COVID-19) pandemic. Despite its quick disappearance worldwide, the strong virulence caused by a few point mutations remains an unsolved problem largely. Along with the other two sublineages, the Delta variant harbors an accumulation of Spike protein mutations, including the previously identified L452R, E484Q, and the newly emerged T478K on its receptor binding domain (RBD). We used molecular dynamics (MD) simulations, in combination with free energy perturbation (FEP) calculations, to examine the effects of two combinative mutation sets, L452R + E484Q and L452R + T478K. Our dynamic trajectories reveal an enhancement in binding affinity between mutated RBD and the common receptor protein angiotensin converting enzyme 2 (ACE2) through a net increase in the buried molecular surface area of the binary complex. This enhanced binding, mediated through Gln493, sets the same stage for all three sublineages due to the presence of L452R mutation. The other mutation component, E484Q or T478K, was found to impact the RBD-ACE2 binding and help the variant to evade several monoclonal antibodies (mAbs) in a distinct manner. Especially for L452R + T478K, synergies between mutations are mediated through a complex residual and water interaction network and further enhance its binding to ACE2. Taking together, this study demonstrates that new variants of SARS-CoV-2 accomplish both "attack" (infection) and "defense" (antibody neutralization escape) with the same "polished sword" (mutated Spike RBD).

Computational modeling of the effect of five mutations on the structure of the ACE2 receptor and their correlation with infectivity and virulence of some emerged variants of SARS-CoV-2 suggests mechanisms of binding affinity dysregulation.

Interactions between the human angiotensin-converting enzyme 2 (ACE2) and the RBD region of the SARS-CoV-2 Spike protein are critical for virus entry into the host cell. The objective of this work was to identify some of the most relevant SARS-CoV-2 Spike variants that emerged during the pandemic and evaluate their binding affinity with human variants of ACE2 since some ACE2 variants can enhance or reduce the affinity of the interaction between the ACE2 and S proteins. However, no information has been sought to extrapolate to different variants of SARS-CoV-2. Therefore, to understand the impact on the affinity of the interaction between ACE2 protein variants and SARS-CoV-2 protein S variants, molecular docking was used in this study to predict the effects of five mutations of ACE2 when they interact with Alpha, Beta, Delta, Omicron variants and a hypothetical variant, which present mutations in the RBD region of the SARS-CoV-2 Spike protein. Our results suggest that these variants could alter the interaction of the Spike and the human ACE2 protein, losing or creating new inter-protein contacts, enhancing viral fitness by improving binding affinity, and leading to an increase in infectivity, virulence, and transmission. This investigation highlighted that the S19P mutation of ACE2 decreases the binding affinity between the ACE2 and Spike proteins in the presence of the Beta variant and the wild-type variant of SARS-CoV-2 isolated in Wuhan-2019. The R115Q mutation of ACE2 lowers the binding affinity of these two proteins in the presence of the Beta and Delta variants. Similarly, the K26R mutation lowers the affinity of the interaction between the ACE2 and Spike proteins in the presence of the Alpha variant. This decrease in binding affinity is probably due to the lack of interaction between some of the key residues of the interaction complex between the ACE2 protein and the RBD region of the SARS-CoV-2 Spike protein. Therefore, ACE2 mutations appear in the presence of these variants, they could suggest an intrinsic resistance to COVID-19 disease. On the other hand, our results suggested that the K26R, M332L, and K341R mutations of ACE2 expressively showed the affinity between the ACE2 and Spike proteins in the Alpha, Beta, and Delta variants. Consequently, these ACE2 mutations in the presence of the Alpha, Beta, and delta variants of SARS-CoV-2 could be more infectious and virulent in human cells compared to the SARS-CoV-2 isolated in Wuhan-2019 and it could have a negative prognosis of the disease. Finally, the Omicron variant in interaction with ACE2 WT, S19P, R115Q, M332L, and K341R mutations of ACE2 showed a significant decrease in binding affinity. This could be consistent that the Omicron variant causes less severe symptoms than previous variants. On the other hand, our results suggested Omicron in the complex with K26R, the binding affinity is increased between ACE2/RBD, which could indicate a negative prognosis of the disease in people with these allelic conditions.

SARS-CoV-2 variant spike and accessory gene mutations alter pathogenesis.

The ongoing COVID-19 pandemic is a major public health crisis. Despite the development and deployment of vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pandemic persists. The continued spread of the virus is largely driven by the emergence of viral variants, which can evade the current vaccines through mutations in the spike protein. Although these differences in spike are important in terms of transmission and vaccine responses, these variants possess mutations in the other parts of their genome that may also affect pathogenesis. Of particular interest to us are the mutations present in the accessory genes, which have been shown to contribute to pathogenesis in the host through interference with innate immune signaling, among other effects on host machinery. To examine the effects of accessory protein mutations and other nonspike mutations on SARS-CoV-2 pathogenesis, we synthesized both viruses possessing deletions in the accessory genes as well as viruses where the WA-1 spike is replaced by each variant spike gene in a SARS-CoV-2/WA-1 infectious clone. We then characterized the in vitro and in vivo replication of these viruses and compared them to both WA-1 and the full variant viruses. Our work has revealed that the accessory proteins contribute to SARS-CoV-2 pathogenesis and the nonspike mutations in variants can contribute to replication of SARS-CoV-2 and pathogenesis in the host. This work suggests that while spike mutations may enhance receptor binding and entry into cells, mutations in accessory proteins may alter clinical disease presentation.

Pandemic-scale phylogenomics reveals the SARS-CoV-2 recombination landscape.

Accurate and timely detection of recombinant lineages is crucial for interpreting genetic variation, reconstructing epidemic spread, identifying selection and variants of interest, and accurately performing phylogenetic analyses1-4. During the SARS-CoV-2 pandemic, genomic data generation has exceeded the capacities of existing analysis platforms, thereby crippling real-time analysis of viral evolution5. Here, we use a new phylogenomic method to search a nearly comprehensive SARS-CoV-2 phylogeny for recombinant lineages. In a 1.6 million sample tree from May 2021, we identify 589 recombination events, which indicate that around 2.7% of sequenced SARS-CoV-2 genomes have detectable recombinant ancestry. Recombination breakpoints are inferred to occur disproportionately in the 3' portion of the genome that contains the spike protein. Our results highlight the need for timely analyses of recombination for pinpointing the emergence of recombinant lineages with the potential to increase transmissibility or virulence of the virus. We anticipate that this approach will empower comprehensive real-time tracking of viral recombination during the SARS-CoV-2 pandemic and beyond.

Insights from a computational analysis of the SARS-CoV-2 Omicron variant: Host-pathogen interaction, pathogenicity, and possible drug therapeutics.

INTRODUCTION: Prominently accountable for the upsurge of COVID-19 cases as the world attempts to recover from the previous two waves, Omicron has further threatened the conventional therapeutic approaches. The lack of extensive research regarding Omicron has raised the need to establish correlations to understand this variant by structural comparisons. Here, we evaluate, correlate, and compare its genomic sequences through an immunoinformatic approach to understand its epidemiological characteristics and responses to existing drugs. METHODS: We reconstructed the phylogenetic tree and compared the mutational spectrum. We analyzed the mutations that occurred in the Omicron variant and correlated how these mutations affect infectivity and pathogenicity. Then, we studied how mutations in the receptor-binding domain affect its interaction with host factors through molecular docking. Finally, we evaluated the drug efficacy against the main protease of the Omicron through molecular docking and validated the docking results with molecular dynamics simulation. RESULTS: Phylogenetic and mutational analysis revealed the Omicron variant is similar to the highly infectious B.1.620 variant, while mutations within the prominent proteins are hypothesized to alter its pathogenicity. Moreover, docking evaluations revealed significant differences in binding affinity with human receptors, angiotensin-converting enzyme 2 and NRP1. Surprisingly, most of the tested drugs were proven to be effective. Nirmatrelvir, 13b, and Lopinavir displayed increased effectiveness against Omicron. CONCLUSION: Omicron variant may be originated from the highly infectious B.1.620 variant, while it was less pathogenic due to the mutations in the prominent proteins. Nirmatrelvir, 13b, and Lopinavir would be the most effective, compared to other promising drugs that were proven effective.

The ctpF Gene Encoding a Calcium P-Type ATPase of the Plasma Membrane Contributes to Full Virulence of Mycobacterium tuberculosis.

Identification of alternative attenuation targets of Mycobacterium tuberculosis (Mtb) is pivotal for designing new candidates for live attenuated anti-tuberculosis (TB) vaccines. In this context, the CtpF P-type ATPase of Mtb is an interesting target; specifically, this plasma membrane enzyme is involved in calcium transporting and response to oxidative stress. We found that a mutant of MtbH37Rv lacking ctpF expression (MtbΔctpF) displayed impaired proliferation in mouse alveolar macrophages (MH-S) during in vitro infection. Further, the levels of tumor necrosis factor and interferon-gamma in MH-S cells infected with MtbΔctpF were similar to those of cells infected with the parental strain, suggesting preservation of the immunogenic capacity. In addition, BALB/c mice infected with Mtb∆ctpF showed median survival times of 84 days, while mice infected with MtbH37Rv survived 59 days, suggesting reduced virulence of the mutant strain. Interestingly, the expression levels of ctpF in a mouse model of latent TB were significantly higher than in a mouse model of progressive TB, indicating that ctpF is involved in Mtb persistence in the dormancy state. Finally, the possibility of complementary mechanisms that counteract deficiencies in Ca2+ transport mediated by P-type ATPases is suggested. Altogether, our results demonstrate that CtpF could be a potential target for Mtb attenuation.

Variants of SARS CoV-2: mutations, transmissibility, virulence, drug resistance, and antibody/vaccine sensitivity.

Severe acute respiratory syndrom coronavirus-2 (SARS CoV-2) is the causative agent of coronavirus disease-19 (Covid-19) which has been designated a worldwide pandemic by the World Health Organization on March 11, 2020. Since that time, the virus has mutated and an assortment of variants have been successful at establishing themselves in the human population. This review article describes the SARS CoV-2 genome, hot spot mutations, variants, and then focuses on the Delta variant, finishing up with an update on the Omicron variant. The genome encompasses 11 open reading frames, one of which encodes the spike or S protein that has been the target for vaccines and some of the drugs because of its role in attachment to the human host cell, as well as antibodies. Mutations in the S protein that are common among several of the variants include D614G that increases transmissibility and viral load and is often associated with P323L on the RNA dependent RNA polymerase. N501Y is a mutation in the receptor binding domain of the S protein that increases binding to the ACE-2 receptor on the human host cells by 10 fold. The discussed variants carry combinations of these and other mutations and are classified by the World Health Organization as variants of concern, variants of interest, and variants under monitoring. All variants are characterized by increased transmissibility (relative to the original SARS CoV-2), which is the reason for their ability to establish themselves. Several but not all variants are more resistant to antiviral drugs and less susceptible to antibodies/vaccines. The Delta variant that dominated the world until November 2021 causes an increased risk for hospitalization and death, but is still very susceptible to the current vaccines. The most recent variant, Omicron, is characterized by increased transmissibility and decreased antibody susceptibility.

The N-Terminal Region of Middle East Respiratory Syndrome Coronavirus Accessory Protein 8b Is Essential for Enhanced Virulence of an Attenuated Murine Coronavirus.

Middle East respiratory syndrome coronavirus (MERS-CoV) is a beta coronavirus that emerged in 2012, causing severe pneumonia and renal failure. MERS-CoV encodes five accessory proteins. Some of them have been shown to interfere with host antiviral immune response. However, the roles of protein 8b in innate immunity and viral virulence was rarely studied. Here, we introduced individual MERS-CoV accessory protein genes into the genome of an attenuated murine coronavirus (Mouse hepatitis virus, MHV), respectively, and found accessory protein 8b could enhance viral replication in vivo and in vitro and increase the lethality of infected mice. RNA-seq analysis revealed that protein 8b could significantly inhibit type I interferon production (IFN-I) and innate immune response in mice infected with MHV expressing protein 8b. We also found that MERS-CoV protein 8b could initiate from multiple internal methionine sites and at least three protein variants were identified. Residues 1-23 of protein 8b was demonstrated to be responsible for increased virulence in vivo. In addition, the inhibitory effect on IFN-I of protein 8b might not contribute to its virulence enhancement as aa1-23 deletion did not affect IFN-I production in vitro and in vivo. Next, we also found that protein 8b was localized to the endoplasmic reticulum (ER)/Golgi membrane in infected cells, which was disrupted by C-terminal region aa 88-112 deletion. This study will provide new insight into the pathogenesis of MERS-CoV infection. IMPORTANCE Multiple coronaviruses (CoV) cause severe respiratory infections and become global public health threats such as SARS-CoV, MERS-CoV, and SARS-CoV-2. Each coronavirus contains different numbers of accessory proteins which show high variability among different CoVs. Accessory proteins are demonstrated to play essential roles in pathogenesis of CoVs. MERS-CoV contains 5 accessory proteins (protein 3, 4a, 4b, 5, 8b), and deletion of all four accessory proteins (protein 3, 4a, 4b, 5), significantly affects MERS-CoV replication and pathogenesis. However, whether ORF8b also regulates MERS-CoV infection is unknown. Here, we constructed mouse hepatitis virus (MHV) recombinant virus expressing MERS-CoV protein 8b and demonstrated protein 8b could significantly enhance the virulence of MHV, which is mediated by N-terminal domain of protein 8b. This study will shed light on the understanding of pathogenesis of MERS-CoV infection.

Increasing the frequency of omicron variant mutations boosts the immune response and may reduce the virus virulence.

The coronavirus has posed a serious threat to the world since its discovery in Wuhan in 2019. Beta, gamma, delta, and the final omicron variants have emerged as a result of several mutations in the virion structure. The Australian Omicron S protein variant contains 37 mutations out of a total of 67 mutations. According to preliminary data from South Africa, Omicron variant infection is not associated with any particular symptoms. The purpose of this research was to determine how changes in the structure of the S protein alter the protein's interaction with the ACE2 receptor. The Omicron variant stimulates the immune response more than the wild strain.

The Genomic Physics of COVID-19 Pathogenesis and Spread.

Coronavirus disease (COVID-19) spreads mainly through close contact of infected persons, but the molecular mechanisms underlying its pathogenesis and transmission remain unknown. Here, we propose a statistical physics model to coalesce all molecular entities into a cohesive network in which the roadmap of how each entity mediates the disease can be characterized. We argue that the process of how a transmitter transforms the virus into a recipient constitutes a triad unit that propagates COVID-19 along reticulate paths. Intrinsically, person-to-person transmissibility may be mediated by how genes interact transversely across transmitter, recipient, and viral genomes. We integrate quantitative genetic theory into hypergraph theory to code the main effects of the three genomes as nodes, pairwise cross-genome epistasis as edges, and high-order cross-genome epistasis as hyperedges in a series of mobile hypergraphs. Charting a genome-wide atlas of horizontally epistatic hypergraphs can facilitate the systematic characterization of the community genetic mechanisms underlying COVID-19 spread. This atlas can typically help design effective containment and mitigation strategies and screen and triage those more susceptible persons and those asymptomatic carriers who are incubation virus transmitters.

Enhanced fitness of SARS-CoV-2 variant of concern Alpha but not Beta.

Emerging variants of concern (VOCs) are driving the COVID-19 pandemic1,2. Experimental assessments of replication and transmission of major VOCs and progenitors are needed to understand the mechanisms of replication and transmission of VOCs3. Here we show that the spike protein (S) from Alpha (also known as B.1.1.7) and Beta (B.1.351) VOCs had a greater affinity towards the human angiotensin-converting enzyme 2 (ACE2) receptor than that of the progenitor variant S(D614G) in vitro. Progenitor variant virus expressing S(D614G) (wt-S614G) and the Alpha variant showed similar replication kinetics in human nasal airway epithelial cultures, whereas the Beta variant was outcompeted by both. In vivo, competition experiments showed a clear fitness advantage of Alpha over wt-S614G in ferrets and two mouse models-the substitutions in S were major drivers of the fitness advantage. In hamsters, which support high viral replication levels, Alpha and wt-S614G showed similar fitness. By contrast, Beta was outcompeted by Alpha and wt-S614G in hamsters and in mice expressing human ACE2. Our study highlights the importance of using multiple models to characterize fitness of VOCs and demonstrates that Alpha is adapted for replication in the upper respiratory tract and shows enhanced transmission in vivo in restrictive models, whereas Beta does not overcome Alpha or wt-S614G in naive animals.

Updated SARS-CoV-2 single nucleotide variants and mortality association.

By analyzing newly collected SARS-CoV-2 genomes and comparing them with our previous study about SARS-CoV-2 single nucleotide variants (SNVs) before June 2020, we found that the SNV clustering had changed remarkably since June 2020. Apart from that the group of SNVs became dominant, which is represented by two nonsynonymous mutations A23403G (S:D614G) and C14408T (ORF1ab:P4715L), a few emerging groups of SNVs were recognized with sharply increased monthly incidence ratios of up to 70% in November 2020. Further investigation revealed sets of SNVs specific to patients' ages and/or gender, or strongly associated with mortality. Our logistic regression model explored features contributing to mortality status, including three critical SNVs, G25088T(S:V1176F), T27484C (ORF7a:L31L), and T25A (upstream of ORF1ab), ages above 40 years old, and the male gender. The protein structure analysis indicated that the emerging subgroups of nonsynonymous SNVs and the mortality-related ones were located on the protein surface area. The clashes in protein structure introduced by these mutations might in turn affect the viral pathogenesis through the alteration of protein conformation, leading to a difference in transmission and virulence. Particularly, we explored the fact that nonsynonymous SNVs tended to occur in intrinsic disordered regions of Spike and ORF1ab to significantly increase hydrophobicity, suggesting a potential role in the change of protein folding related to immune evasion.

Enhanced fusogenicity and pathogenicity of SARS-CoV-2 Delta P681R mutation.

During the current coronavirus disease 2019 (COVID-19) pandemic, a variety of mutations have accumulated in the viral genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and, at the time of writing, four variants of concern are considered to be potentially hazardous to human society1. The recently emerged B.1.617.2/Delta variant of concern is closely associated with the COVID-19 surge that occurred in India in the spring of 2021 (ref. 2). However, the virological properties of B.1.617.2/Delta remain unclear. Here we show that the B.1.617.2/Delta variant is highly fusogenic and notably more pathogenic than prototypic SARS-CoV-2 in infected hamsters. The P681R mutation in the spike protein, which is highly conserved in this lineage, facilitates cleavage of the spike protein and enhances viral fusogenicity. Moreover, we demonstrate that the P681R-bearing virus exhibits higher pathogenicity compared with its parental virus. Our data suggest that the P681R mutation is a hallmark of the virological phenotype of the B.1.617.2/Delta variant and is associated with enhanced pathogenicity.

Identification of genotypic variants and its proteomic mutations of Brazilian SARS-CoV-2 isolates.

The second wave of COVID-19 caused by severe acute respiratory syndrome virus (SARS-CoV-2) is rapidly spreading over the world. Mechanisms behind the flee from current antivirals are still unclear due to the continuous occurrence of SARS-CoV-2 genetic variants. Brazil is the world's second-most COVID-19 affected country. In the present study, we identified the genomic and proteomic variants of Brazilian SARS-CoV-2 isolates. We identified 16 different genotypic variants were found among the 27 isolates. The genotypes of three isolates such as Bra/1236/2021 (G15), Bra/MASP2C844R2/2020 (G11), and Bra/RJ-DCVN5/2020 (G9) have a unique mutant in NSP4 (S184N), 2'O-Mutase (R216N), membrane protein (A2V) and Envelope protein (V5A). A mutation in RdRp of SARS-CoV-2, particularly the change of Pro-to Leu-at 323 resulted in the stabilization of the structure in BRA/CD1739-P4/2020. NSP4, NSP5 protein mutants are more virulent in genotype 15 and 16. A fast protein folding rate changes the structural stability and leads to escape for current antivirals. Thus, our findings help researchers to develop the best potent antivirals based on the new mutant of Brazilian isolates.